Transport coefficients from the 2PI effective action 
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We show that the lowest nontrivial truncation of the two-particle irreducible (2PI) effective action 
correctly determines transport coefficients in a weak coupling or 1/N expansion at leading (loga- 
rithmic) order in several relativistic field theories. In particular, we consider a single real scalar 
field with cubic and quartic interactions in the loop expansion, the O(N) model in the 2PI-1/AT ex- 
pansion, and QED with a single and many fermion fields. Therefore, these truncations will provide 
a correct description, to leading (logarithmic) order, of the long time behavior of these systems, 
i.e. the approach to equilibrium. This supports the promising results obtained for the dynamics of 
quantum fields out of equilibrium using 2PI effective action techniques. 
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Introduction. Recent developments in heavy-ion col- 
lisions and cosmology have spurred the theoretical un- 
derstanding of the dynamics of quantum fields out of 
equilibrium. In particular, the thermalization of quan- 
tum fields is a subject of both fundamental and practi- 
cal relevance. For quantum fields far from equilibrium, 
promising results have been obtained from a systematic 
use of the 2PI effective action 0, formulated along the 
Schwingcr-Keldysh contour. While the basic formulation 
of this approach is well-known (see Refs. for re- 
cent applications in equilibrium) , the recent progress has 
been the numerical solution of the resulting evolution 
equations for the (one- and) two-point functions with- 
out any further approximation. This allows one to go 
far from equilibrium and describe, e.g., the emergence of 
quasi-particles in a completely self-consistent way. This 
program has been carried out for a single scalar field with 
quartic self-interactions using a three-loop expansion in 
1 + 1 dimensions 0, for the O(N) model using the 2PI- 
1/N expansion in 1 + 1 00 and 3+1 dimensions 0, 
and recently also for a chirally invariant Yukawa model 
in 3 + 1 dimensions . The (mostly numerical) results 
obtained so far suggest that the lowest nontrivial trun- 
cation beyond the mean-field approximation is necessary 
and sufficient to describe in one formalism both the dy- 
namics far from equilibrium as well as the subsequent 
equilibration. 

A necessary requirement for any method to success- 
fully describe nonequilibrium field dynamics is that it en- 
compasses the correct long-time behavior. Close to equi- 
librium, the evolution of the system on long time- and 
length scales is characterized by transport coefficients, 
which have been computed at leading order in a weak 
coupling or 1 /N expansion [lCl IllL [l2j • In order to as- 
sess the validity of truncations of the 2PI effective action, 
it is therefore crucial that transport coefficients obtained 
within the 2PI formalism agree, in the weak coupling 
limit, with those results. 

In this paper, we show how the calculation of trans- 
port coefficients is organized in the framework of the 2PI 



effective action. We then consider a variety of theories 
and show, by comparing with results obtained previously, 
that truncations currently used in far-from-equilibrium 
applications include in the weak coupling or large N 
limit the appropriate diagrams to yield the correct result 
for transport coefficients to leading (logarithmic) order. 
This result provides strong support for the applicability 
of truncations of the 2PI effective action to describe the 
equilibration of quantum fields. 

2PI effective action. The Kubo formula relates trans- 
port coefficients to the expectation value of appropriate 
composite operators in thermal equilibrium (to be pre- 
cise, to the imaginary part of retarded correlators at zero 
momentum and vanishing frequency). The correlator we 
will be interested in here is therefore of the form 



(0(x)0(y))-(0(x))(0(y)), 



(1) 



with O an operator bilinear in the fundamental fields. For 
example, for the shear viscosity 0{x) = Uij(x), where for 



a real scalar field 7T; 



Or 



^5ijdk4>dk4>, and for the 
electrical conductivity 0(x) = j z (x), with j l = $-f l 4>. 

We first demonstrate that the 2PI effective action gen- 
erates precisely correlators of the form (TIJ . We consider 
the case of a real field <fi, coupled to a bilinear source K 
(the extension to fermions is straightforward). The path 
integral is 



(2) 



with S the classical action. (We use a condensed no- 
tation where latin indices summarize space-time as well 
as internal indices, and integration and summation over 
repeated indices is understood.) For the application 
we discuss here it is not necessary to couple a source 
to <j) itself and we assume throughout that (cf> 1 ) van- 
ishes. The 2PI effective action is defined as the Leg- 
endre transform of W, T[G] = W[K] - \G i] K tj , with 
8W[K]/5Kij = |G«, and 8T[G)/8G^ = -\Kij. The 
two-point function G y = (Tccffft) is the (connected) 
time-ordered two-point function along the contour C in 
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FIG. 1: Expectation value of bilinear operators (black dots) 
from the 2PI effective action. The shaded square denotes the 
connected 4-point function (with external legs amputated). 
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FIG. 2: Integral equation for the 4-point function, derived 
from the 2PI effective action. 



the complex-time plane. Our discussion will be general 
and we do not need to specify the contour here; at any 
moment one may specialize to the Matsubara contour or 
the Schwinger-Keldysh contour. 

In order to obtain a correlator of the form (JIJ , we dif- 
ferentiate W twice with respect to K: 



8 2 W[K] 
8K tj 5K kl 



- [{Tc^V^V) - G ij G kl ] 



4 L c 



G G 



(3) 



where G l J' ,kl is the usual connected 4-point function 
(semicolons separate indices with a different origin). We 
note that W does not generate connected Green func- 
tions, instead it generates a 4-point function which, after 
identifying i,j with x and k, I with y, is precisely of the 
form (see Fig.^). To proceed further we remove the 
external legs, Gf' M = G ijf G"' G kk ' G 11 ' vf], . k , v , and con- 
centrate on the 4-point vertex function r' 4 ) . Note that in 
a theory with cubic interactions is not 1-particle irre- 
ducible. The vertex function obeys an integral equation 
that can be obtained using standard functional relations 
Q. It reads (see Fig. 01 



r< 4 > - A 

1 ij;kl — ix irM 



(4) 



To apply the Kubo formula, we specialize to a system 
in thermal equilibrium that is invariant under spacetime 
translations. In momentum space, Eq. Q then reads 

T^(p,k)=A(p,k) + ^J A(p, q)G 2 (q)T W (q,k). (7) 

The importance of using Eq. Q in the renormalization of 
the gap equation has been emphasized recently 1 1 ll] . 
Eq. Q is valid both in the imaginary-time as well as in 
the real-time formalism, where the 4-point function and 
kernel have a more complicated tensor structure |l5|. 

We note that the 4-point function appears in the inte- 
gral equation with a particular momentum configuration: 
the momentum p (k) enters and leaves on the left (right) 
and the two intermediate propagators carry the same 
momentum q. This configuration suffers therefore from 
pinching poles: when the loop momentum q is nearly on- 
shell, the product of propagators is potentially very large 
and all terms in the ladder series may be equally im- 
portant. This situation is precisely the one that appears 
in the diagrammatic evaluation of transport coefficients 
[irj|: computing transport coefficients diagrammatically 
amounts to studying an integral equation of the type JJJ, 
specialized to the case that the external momenta p and k 
as well as the internal momentum q are on- shell liana. 

Loop expansion. As a first example, we consider a real 
scalar field with cubic and quartic interactions, 



S = 



-m 2 (f> 2 



fJ 

3!< 



4! 



(8) 



In the loop expansion for r2[G], we include terms up to 
three loops (see Fig. I2J . The self energy £ is obtained 
by cutting the diagrams once. Cutting the result once 
more yields the kernel A in the integral equation for the 
4-point function (see Fig.^J. 

We now show that this truncation includes the physics 
relevant for the leading-order result of the shear viscos- 
ity in the weak-coupling limit, A ~ (g/m) 2 < 1. In 
order to do this, it is sufficient to show that Eq. (JTJ 



where the 4-point kernel follows from 



A 



ij ; kl 



S^iG] 



2i 



ST 2 [G] 



SG kl ' lJ " 5Gv 
when the effective action is written as 111 



(5) 



T[G] = ^TrlnG- 1 + ^TrGo \G - G ) +T 2 [G]. (6) 

Here, iG$ is the free inverse propagator and T 2 [G] is the 
sum of all 2-particle irreducible (2PI) diagrams with no 
external legs and exact propagators on the internal lines. 
We emphasize that Eq. Q is exact and applies out of 
equilibrium. The first derivative of T[G] determines the 
gap equation G _1 = Gq 1 — E in the absence of sources. 




FIG. 3: Contributions to the 2PI effective action in the loop 
expansion up to three loops in a theory with cubic and quartic 
interactions. 
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FIG. 4: Zero- and one-loop skeleton kernel from the contri- 
butions to the 2PI effective action shown in Fig. |5] 



FIG. 5: One-loop perturbative kernel from the 2PI effective 
action at two-loop order in a theory with cubic interactions. 
The dashed line indicates how to cut the diagrams to make 
the connection with the scattering kernel in kinetic theory. 



includes all diagrams in the appropriate kinematic con- 
figuration which are known to contribute at leading or- 
der 0] ■ We therefore specialize to on-shell momenta p, 
k, and q, where the leading pinching-pole limit arises. 
One must then proceed to a careful analysis of all the 
possible perturbative contributions to the integral equa- 
tion from the kernel A(p, k). For example, the single rung 
with a bare propagator does not contribute straightfor- 
wardly due to kinematics, but it does when either it is 
iterated in the integral equation to get a one-loop kernel 
or a one-loop self-energy correction is included, as de- 
picted in Fig . (for detailed power-counting arguments, 
see Ref. [lOj ) . In general, Eq. J7J will contain contri- 
butions that are subleading in the weak-coupling limit. 
For instance, the four final rungs in Fig. ^ can be seen 
as vertex corrections to the single rung and contribute 
at subleading order only. When the result of the power- 
counting analysis of Ref. is carried over to this case, 
one finds that the leading-order contribution to the ker- 
nel can be written as an integral over an effective scat- 
tering kernel \X + g 2 [Gr(s) + Gn(t) + Gr(u)]\ , where 
we write s, t, u to indicate the contributions from the 
three scattering channels and Gr denotes the retarded 
Green's function. This kernel is the square of the sum of 
all 2-to-2 processes in a theory with cubic and quartic in- 
teractions, which establishes the connection with kinetic 
theory 10]. It is instructive to see how this result is put 
together. Consider for a moment only the two-loop dia- 
gram with cubic interactions (see Fig.EJ. The scattering 
kernel that arises from this diagram reads explicitly (see 
e.g- 0): g A [\GR(s)\ 2 + \G R (t)\ 2 + \G R (u)\%i.e.,theswn 
of the squares of the matrix elements, however, without 
interference terms. Indeed, it is easy to see that the in- 
terference terms originate from the 3-loop diagrams. 

We conclude therefore that the 3-loop approximation 
of the 2PI effective action is necessary and sufficient to 
yield the leading-order result for the shear viscosity in a 
weakly-coupled scalar theory [lfij ]. 

O(N) model. We now consider a real scalar TV- 
component quantum field </> a (a = 1, ...,N) with a 
classical 0(-/V)-invariant action and an interaction term 
(X/4:\N)((j) a (j) a ) 2 . Instead of the loop expansion, we 
will use the 2PI-l/iV expansion to next-to-leading or- 
der (NLO), which is discussed in detail in Refs. 0, @- 
The 2PI part of the effective action can be written as 



FIG. 6: Contributions to the 2PI effective action in the O(N) 
model in the 2PI-l/iV expansion at LO and NLO. Only the 
first few diagrams at NLO are shown. 
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FIG. 7: Integral equation for the 4-point function in the 2PI- 
1/N expansion of the O(N) model at NLO. 
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FIG. 8: Integral equation for the auxiliary correlator D. 



T 2 [G] = r%°[G] + T™ LO [G] + ..., with (see Fig.©) 
r^°[G] = J G aa {x,x)G bb {x,x), 



r^ LO [G] 



-TrlnB(G), 



(9) 
(10) 



where B(x,y;G) = 6 c {x-y) + (iX/6N)G ab (x ) y)G ab (x ) y) 
sums bubbles (which can be seen by re-expanding the 
logarithm) . In this case the kernel reads (see Fig. EJ) 



, NLO 



(11) 



Kb,cd(P,k) = -[8acSbd + S ad d bc \D(p-k) 



+2 / G ab (p - q)D 2 (q)G cd (k - q), (12) 

Jq 

where we used the auxiliary correlator D = 
(«A/3A r )B _1 to sum the chain of bubbles (see 
Fig. EJ: 



D{p) = 



Ii{p)D{p)] 1 with 



(iX/'SN) [1 
IL(p) = -lf q G ab (q)G ab (p + q) [|. 

The advantage of the large N expansion employed here 
is that it allows for a computation of the shear viscosity 
to first nontrivial order in the 1/N expansion without a 
restriction to small A. This is especially pressing for ap- 
plications of the 0(4) model to QCD phenomenology in 
which the coupling constant A has to be taken large. As 
far as we know, such a calculation of the shear viscosity 
has not yet been performed. The integral equation de- 
rived here provides a convenient starting point for that 
and work in this direction is currently underway [l7| . 

QED. As a last example we consider QED. Although 
for gauge theories there are a number of issues related 
to the gau ge s ymmetry (dependence on the gauge-fixing 
parameter |l8j. Ward identities 0]) which have not been 
resolved completely, we think it is nevertheless worth an- 
alyzing to what accuracy transport coefficients can be 
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a) b) c) 

FIG. 9: Contributions to the 2PI effective action in QED in 
the loop expansion with (a) 2 and (b) 3 loops or in the 2PI- 
1/N expansion at (a) NLO and (b,c) NNLO. 
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FIG. 10: Integral equations for the 4-point functions in QED 
at 2-loop order in the loop or at NLO in the 1/N expansion. 



expected to be computed within a specific truncation of 
the 2PI effective action. 

In the loop expansion (see Fig.^J, we restrict the dis- 
cussion to the 2-loop approximation. The corresponding 
integral equations are shown in Fig. 1101 In order to an- 
alyze what processes are included in the weak-coupling 
limit, we follow the strategy of the scalar theory and 
rewrite the integral equation such that the external mo- 
menta are on-shell. The possible one-loop contributions 
to the kernel are shown in Fig. ^2 This kernel contains 
the sum of the squares of all 2-to-2 processes, but not the 
interference terms. Indeed, it is the 3-loop diagram that 
is responsible for interference. Using the results from 
the diagrammatic analysis [2fJ, we find that, in order to 
determine the shear viscosity and the electrical conduc- 
tivity to leading-log order, only the first two rungs in 
Fig. Illl need to be considered, with the momentum flow- 
ing through the rung being soft and below the lightcone. 
We conclude that the integral equations in the 2-loop 
approximation sum the necessary diagrams [2(1 \2l\ to 
obtain the shear viscosity and electrical conductivity in 
QED to leading-logarithmic order. 

Finally, we consider QED with N identical fermions 
and a rescaled coupling e 2 — > e 2 /N in the large N limit 
(see Fig. We consider only the NLO contribution. 




FIG. 11: One-loop perturbative kernel in the 2-loop approxi- 
mation in QED. In the 1/N expansion, the first and the fourth 
diagrams dominate. 



In this case the kernel is dominated by the first and the 
fourth diagrams in Fig. 1111 the other contributions are 
suppressed by 1/N. The photon propagator still resums 
fcrmion bubbles and is in fact similar to the auxiliary cor- 
relator in the O(N) model. Cutting the first and fourth 
diagrams as in Fig. \5\ shows that they correspond to 
Coulomb scattering in all three channels, which are in- 
deed 0| the dominant processes in a kinetic theory with 
many fermions. We conclude that the 2PI-1/7V expan- 
sion in QED with N fermions at NLO would yield the 
correct leading-order result for the shear viscosity and 
electrical conductivity. 

In summary, we have shown how the calculation of 
transport coefficients is organized in the framework of 
the 2PI effective action. For a variety of models, we 
have discussed the first nontrivial truncations in a weak 
coupling or 1 /N expansion and found that these trunca- 
tions yield ladder integral equations with the particular 
kinematic configuration and rungs appropriate to obtain 
transport coefficients at leading (logarithmic) order in 
the weak coupling or 1/N expansion. This formulation 
offers therefore a systematic starting point for the deriva- 
tion of these integral equations. Methods for the actual 
solution of the integral equations (a topic not discussed 
here) can be found in the literature. 

In the wider context of nonequilibrium quantum field 
theory, our findings provide theoretical support for the 
successful description of quantum fields out of equilib- 
rium using truncations of the 2PI effective action, which 
are currently actively under investigation. We would like 
to emphasize that any scheme which aspires to properly 
describe the nonequilibrium evolution and ensuing ther- 
malization of a system has to render transport coefficients 
correctly. 
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